Weak localization of disordered quasiparticles in the mixed superconducting state 
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Starting from a random matrix model, we construct the low-energy effective field theory for the 
noninteracting gas of quasiparticles of a disordered superconductor in the mixed state. The theory 
is a nonlinear a model, with the order parameter field being a supermatrix whose form is determined 
solely on symmetry grounds. The weak localization correction to the field-axis thermal conductivity 
is computed for a dilute array of s-wave vortices near the lower critical field H&- We propose that 
weak localization effects, cut off at low temperatures by the Zeeman splitting, are responsible for 
the field dependence of the thermal conductivity seen in recent high-T c experiments by Aubin et al. 
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I. INTRODUCTION 

The long wave length physics of phases of matter with 
spontaneously broken symmetries is commonly described 
by an effective field theory for the relevant order param- 
eter. For the problem of localization and transport in 
disordered metallic systems at low temperature, the ap- 
propriate "order parameter" is known Q to be a super- 
matrix (or a matrix of dimension zero if the replica trick 
is used), conventionally denoted by Q. Three univer- 
sality classes, differing by their behavior under time re- 
versal and spin rotations, are widely known H. They 
are labeled by an index j3 — 1,2,4, and are traditionally 
referred to as the classes with orthogonal, unitary, and 
symplectic symmetry. We denote them by Al, A, and All 
for short |J. In each case, the field theory for Q belongs 
to the general family of nonlinear a models. The field Q 
contains the Goldstone modes of a hidden symmetry 
connecting retarded and advanced single-electron Green 
functions, which is broken by a nonzero density of states. 
At tree level one recovers the classical diffusion approxi- 
mation, which neglects quantum interference corrections 
due to electron paths with loops. Anharmonic terms in 
the field theory represent "interactions" between the dif- 
fusion modes, giving rise to so-called weak localization 
corrections to diffusion. For the classes Al and A in di- 
mension d < 2, these interactions become strong at large 
distance scales and thus cause localization of all states, 
regardless of the strength of the disorder. 

In recent years it was found that the (3 = 1,2,4 classifi- 
cation is not exhaustive: for systems with symmetries of 
the particle-hole (ph) type, the invariance group of the 
order parameter field Q becomes enlarged in the vicin- 
ity of the ph-symmetric point. One instance of such 
symmetry enhancement are Dirac fermions in a random 
gauge field J|,||, another are disordered quasiparticles 
that exchange charge (but no energy) with a supercon- 
ducting condensate || . Such systems exhibit novel spec- 
tral statistics and transport properties. Seven symmetry- 
enhanced universality classes have been identified, and 
the corresponding order parameters Q were constructed 



from their random matrix limit in Q . The main message 
is that Q lives on a symmetric space - in precise techni- 
cal language: on a Riemannian symmetric superspace - 
in all cases. The nonlinear a models defined over such 
spaces are known jg] to be attractive under the flow of 
the renormalization group. Therefore the order param- 
eter of a given disordered single-particle system, and in 
most cases its low-energy effective field theory as well, 
can be inferred quite simply by investigating the ergodic 
(or random matrix) limit. 

The present report focusses on class C, which emerges 
for noninteracting low-energy quasiparticles in a mag- 
netic field and in contact with a spin-singlet supercon- 
ductor. The denning condition || is that the quasipar- 
ticle Hamiltonian be invariant under SU(2) rotations of 
the electron spin, whereas time reversal invariance has 
to be broken. Since a superconductor screens magnetic 
fields, this universality class can only be realized in an 
inhomogeneous superconducting state - unless time re- 
versal invariance is broken spontaneously. In the very 
recent literature the following realizations have appeared 
Q: i) a metallic quantum dot in the form of a chaotic 
billiard, subject to a magnetic flux and bordering on a 
superconductor pfj| , pT[ ; ii) quasiparticles in the core of 
an isolated vortex in a disordered s-wave superconductor 
p2[ ; and iii) a (quantum) disordered version of a d x 2_ y 2 
superconductor with orbital coupling to a magnetic field 
|fl3f . The hallmark of class C is that, in contrast with the 
metallic class A, the weak localization correction does not 
vanish Jl^l , in spite of the presence of a magnetic field. 
The persistence of weak localization in a field is caused 
by nonstandard modes of quantum interference that ap- 
pear when impurity and Andreev scattering are simulta- 
neously present. In a semiclassical picture, the effect can 
be understood as being due || to quasiparticle paths in 
which a loop is circled twice, with the charge states dur- 
ing the first and second looping being exactly opposite 
to each other. 

To identify the order parameter field Q and its low- 
energy effective theory for class C, one may proceed in 
several ways. The direct method, due to pa and worked 
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out in detail for isolated vortices of an s-wave super- 
conductor in Ref. is to start from the BCS mean 
field Hamiltonian for the quasiparticles, set up a super- 
symmetric generating functional for the Gorkov Green 
function, introduce a composite field Q to decouple the 
4-vertices produced by averaging over the disorder, inte- 
grate out the quasiparticle fields, solve two saddle point 
equations for Q in sequence (the second of which turns 
out to coincide with the Usadel equation [jl6| ) , and finally 
expand in gradients of Q to obtain the low-energy effec- 
tive theory. The field theory so obtained is a nonlinear 
a model, with Q taking values in a Riemannian symmet- 
ric superspace of type DIIIICI, in agreement with the 
random matrix analysis of M. Its coupling constant has 
the universal meaning of a conductivity for the conserved 
probability (or energy) current transported by the quasi- 
particles. Because quasiparticles also carry spin, the cou- 
pling constant may be reinterpreted jl3) as a spin conduc- 
tivity in the present context. (The latter interpretation 
fails for systems with spin-orbit scattering or magnetic 
impurities, where spin is not conserved.) 

Given the proper identification of the order parame- 
ter field Q, a few qualitative conclusions are immediate. 
According to the renormalization theory of nonlinear a 
models B , the sign of the one-loop renormalization group 
beta function in two space dimensions is completely de- 
termined by the sign of the (Ricci) curvature tensor rel- 
ative to the metric tensor. Since the curvature of the 
Riemannian symmetric superspace of type Dili | CI is 
positive |l7]] , the weakly coupled two-dimensional theory 
renormalizes by logarithmic corrections towards strong 
coupling (i.e. strong disorder), which ultimately leads 
to localization of all quasiparticle states at T = 0. This 
localized phase was called a "spin insulator" in |L3|. In 
dimension d — 1 the same corrections are present, but 
with a linear dependence on the cutoff length. In d = 3 
the theory supports a delocalization transition to a phase 
of extended states, the "spin metal" The addition 

of random classical Heisenberg impurity spins (at sub- 
critical concentration, so as to maintain superconductiv- 
ity) causes crossover from class C to class D ||. with 
the nonlinear a model changing to type CI|Z?III |Q. In 
the process, the sign of the symmetric space curvature 
gets reversed, whence weak localization turns into weak 
antalocalization, making it possible for extended states 
to exist already in two dimensions. 

Our goal here is to extend the treatment of |l^] and 
illustrate some of the above general facts at the thermal 
transport of the class C quasiparticles of a disordered 
s- or d-wave superconductor in the mixed state [T^ ] . We 
assume the presence of (nondescript) nonmagnetic impu- 
rities, which disorder the vortex array and cause elastic 
scattering of the quasiparticles. To tackle this problem 
we will use a coarse grained or random matrix type of 
approach, placing the emphasis on symmetry considera- 
tions. 



II. EFFECTIVE FIELD THEORY FROM AN 
iV-ORBITAL MODEL 

We begin our treatment by partitioning the supercon- 
ductor into cells of equal size, with each cell containing 
one vortex segment with a length of the order of the elas- 
tic mean free path I. Within each cell we introduce (in 
the spirit of the real-space renormalization group) a basis 
of N quasiparticle wavefunctions that comprise the rele- 
vant low-energy configurations. The matrix of the Hamil- 
tonian in such a basis assumes a sparse block structure, 
with one block on the diagonal for each cell, and with 
off-diagonal blocks that couple neighboring cells. If i la- 
bels the cells and a — 1,...,N the orbitals inside a cell, 
the "coarse grained" Hamiltonian is of the form 

H = \ h ™,iMa] C m + C lal C 3 bl) 

+Aia,j6(c| aT ct H -ct ai ct bT )/2 + h.C.) , 

where the sum over i,j is restricted to i = j and pairs of 
neighboring cells. The spin-singlet nature (fj. — J.j) of 
the coupling to the pairing field is dictated by conserva- 
tion of spin. Fermi statistics then requires the complex 
matrix A to be symmetric: lS.ia.jb = ^-jb,ia p9| - If we 
temporarily suppress the cell and orbital indices, H can 
be written in the schematic form H = Tr(7Ycc) + const, 
where 

-CD. e-M,, «-(*; «). 

The symmetries of the Hamiltonian matrix TL are sum- 
marized by the equation TL — —CTL T C~ 1 , with C being 
the symplectic unit C — io~2 ®\. Note that when the Zee- 
man energy H z = M s Eia( c Lt c iaT _ 4ai c ia.i)/ 2 is taken 
into account, the SU(2) spin rotation invariance of H is 
broken down to a U(l) symmetry. 

Disorder in the microscopic Hamiltonian gives rise to 
randomness in TL. Because the universal properties at 
long wave lengths are insensitive to the microscopic de- 
tails, we have considerable freedom in choosing the ran- 
dom Hamiltonian TL. The simplest choice is an iV-orbital 
model with locally gauge-invariant disorder of the type 
invented by Wegner [^0) for the purpose of describing the 
universal physics of the Anderson localization transition 
for (3 = 1,2, 4. The crucial new feature in the present case 
is the relation TL = — C7i T C _1 , which is invariant under 
symplectic transformations TL i— > STLS^ 1 , S T CS = C. 
We therefore adopt a model with local Sp(2A^) gauge in- 
variance: the elements of the matrix TL are taken to be 
Gaussian distributed uncorrelated random variables with 
zero mean, (TL) = 0, and second moments specified by 

(TrATLij TrBTL kl ) - ^ (S^TrAB - 5$ TiACB^C ~ 1 ) 
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where Sfj = SikSji, and Wij is a rapidly decreasing func- 
tion of the distance between the cells i and j. Aside from 
respecting the symmetries and locality of the Hamilto- 
nian, this choice has the virtue of maximizing the in- 
formation entropy. The main benefit from using such a 
maximum entropy model is that the introduction of the 
supermatrix Q, usually a tricky step that requires some 
expertise, becomes straightforward as we now proceed to 
show. 

We replace the operators c,c by classical fields ip,ip: 
CiaaCjbp — > J2<7 i>iaa,<Tip<r,jbp and integrate bilinears in 
ip, ip against expiTr(7i — E)ipip in the usual way to gen- 
erate the Gorkov Green function at energy E. The intro- 
duction of a bosonic partner (a = B) for each fermionic 
field (<r = F) serves to cancel vacuum graphs by the 
mechanism of supersymmetry. There is one complica- 
tion, however: the matrix iptp does not share the symplec- 
tic symmetry of the Hamiltonian. To remedy this mis- 
match, we introduce an extra quantum number ( "pseudo 
charge") c = ±1, so that the quasiparticle fields expand 
to tensors ipiaa,<7c and ip ac ,iaa M- On imposing the con- 
ditions ijj = Cip T "f^ 1 and ip — — 7?/> T C~ 1 , where 7 is a 
real orthogonal matrix that will be specified shortly, we 
have the symmetry ipip = — C(^-0) T C _1 as desired. 

Since the order parameter Q is a local field, its nature 
can be uncovered by looking at the Hamiltonian trun- 
cated to a single cell. With this truncation temporarily 
in force, we introduce Q as follows: 



dHexp(-NTrH 2 /2wa + iTiH^ip) 
dQexp(-NSTrQ 2 /2w + iSTrQVVO . 



The equality is verified by using the cyclic invariance of 
the (super)trace: Tr(>V>) 2 = STr(V^) 2 . The Hubbard- 
Stratonovitch field Q is a 4 x 4 supermatrix which, by its 
coupling to iftp, inherits the symmetry 



Q = -iQ T r 



The constraints relating ip and ip to one another are com- 
patible only if 7 2 equals the superparity matrix (+1 on 
bosons, —1 on fcrmions). To meet this condition we put 
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E, 



BB ® CTi 



E 



FF 



ia 2 



0"i 






icr 2 



Relation ([!]) is the defining equation of an orthosymplec- 
tic Lie algebra and is invariant under Q t— > TQT~ X with 



7(T" 



-1\T 
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T e OSp(2|2). In the general case, where 



Green functions at n different energies are to be averaged, 
Q acquires matrix dimension 4u x 4n, and the symmetry 
group gets enlarged to OSp(2n|2n) = G. 

Returning now to the full lattice problem, introducing 
Qi for every cell i and integrating over the quasiparticle 



fields ip,ip we arrive at the following action functional: 
S/N = w- 1 lJ STrQ l Q J /2 + ^ STr In(Q ( - uo ® £3) 



where {^^aca'd = 5 aa '(a^) cc i and w is a diagonal ma- 
trix containing the energies at which the quasiparticle 
Green functions are to be evaluated. Variation of S yields 
the saddle point equation w~ 1 ijQj = (^£3 — Qi) -1 , 
whose physical solution (dictated by causality of the 
Green function) at lo — and homogeneous in space 



is Q° — ivT,3 with v 



Low-energy fluc- 



tuations result from setting Qi = TiQ°T~ 1 and taking 
Ti E G to vary slowly with the position of the cell i. By 
expanding in gradients, the low-energy effective action 
for such configurations at lo = is easily seen to be a 
nonlinear a model, 
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j d 3 xSTr(z? ± (V ± g) 2 +D||(V||Q) 2 ) , (2) 



where we have switched to continuous coordinates x±_ = 
(x,y) and xy = z. The parameter v is the density of 
states of the superconductor, and Dn , D± are the field- 
axis and transverse effective diffusion constants of the 
quasiparticle gas. We are using units K = 1. At finite u>, 
the field theory action is perturbed by a term 



S u = — — I d 3 x STrwSaQ 



S — Sq + S u 



(3) 



We have rescaled the field to Q = TE 3 T _1 . Since this 
expression for Q is invariant under T — > Tk for k = 
Y>j,kYj3 G GL(n|n) = K, the supermatrix Q lives on a 
coset space G/K. If we parametrize Q by 



exp 



X 
X 



positivity of So or equivalently, stability of the functional 

t 



(V\ integral, requires Abb = +A BB and Aff 



FF' 

In invariant mathematical language, this means Qbb £ 
S0*(2n)/U(n) |l) and Q FF G Sp(2n)/U(n), which are 
symmetric spaces of type Dili and CI - hence the name 
Dili I CI for the present nonlinear a model. The same ef- 
fective theory (restricted to the FF sector due to the use 
of fermionic replicas) was obtained in Ref. |13|, based on 
a quasiparticle Hamiltonian for a dirty d x i_ y i supercon- 
ductor with orbital coupling to a magnetic field. This is 
no surprise, as that system belongs to symmetry class C 
and the order 'parameter field Q and its low-energy effec- 
tive theory are determined solely by symmetry. (Inciden- 
tally, the classification scheme of || assigns the quasipar- 
ticles of the d x 2_ y 2 superconductor in zero field to class 
CI. According to Q, the corresponding symmetric su- 
perspace is D|C, also in agreement with the findings of 
Ref. @.) 
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In order to break parity and account for the Hall an- 
gle, one would need to add to the Lagrangian a topo- 
logical density proportional to e kl STrQdkQdiQ, which 
is closely related to Pruisken's term p3] well known 
from the theory of the integer quantum Hall effect. In 
two dimensions this term integrates to a winding number 
and is nontrivial, since the fundamental group of U(n) is 
ni(U(n)) = Z and there exists the topological identity 
ni(U(n) = n 2 (Sp(2n)/U(n)). However, such a topolog- 
ical term does not affect the results for the longitudinal 
spin and thermal conductivities presented below and will 
therefore be omitted. 

The maximum entropy derivation presented here does 
not supply microscopic expressions for the couplings vD\\ 
and uD± . (We can express them in terms of the random 
matrix parameters Wij, but this is neither illuminating 
nor useful.) These parameters can either be calculated 
from (quasi)classical transport theory |23| or, better yet, 
taken from experiment. In the latter case we extract the 
(bare) coupling constants from experiments conducted at 
temperatures high enough so that the transport is clas- 
sical, and then use the field theory (||,||) to predict the 
quantum corrections that emerge at lower temperatures. 



III. WEAK LOCALIZATION IN CLASS C 

The field theory (|^) does not apply to charge trans- 
port, as the condensate carries charge and quasiparticle 
charge is not a constant of motion. The energy, how- 
ever, and for class C also the spin of a quasiparticle are 
conserved, which allows to probe for quasiparticle trans- 
port and localization by measuring the thermal and spin 
transport. To obtain the relevant transport coefficients 
we start from the bilocal conductivity tensor 



tu 



(x,x';£) =^«WGl,(x,x' i £)«WG^( X ' ) xi£) ) 



which describes the nonlocal linear response of the spin 
current to a perturbation due to the Zeeman coupling 
with an applied field. The quantities G R and G A are 
the retarded and advanced Gorkov Green functions and 
vi a = (i(o'3)aa(di — di) — 2eA;)/2m is the Z-component 
of the velocity operator. We use the relation H = 
—CTf^C^ 1 to express G A at energy E by G R at energy 
—E |p^| . Disorder averaging and the mapping on the 
nonlinear a model with n = 2 then turn the tensor tu into 
a correlation function of the conserved OSp(4|4) Noether 
current J = (QVQ) B11322 of the field theory: 



(r«(x,x')) = (™A) 2 (JKx)^(x')) 



(4) 



The second superscript in the expression for the Noether 
current refers to the pseudo charge, while the third dis- 
tinguishes between the two Green functions. The sym- 
metry breaking perturbation due to the quasiparticle en- 



ergy E is incorporated into the formalism by setting 
uj = diag(£+, -E~), where E ± = E ± iO. 

Next, let an = ^tt)" 1 / (r H (x, x')) d 3 x' be the local 
"spin" conductivity for quasiparticles with fixed spin up 
or down. To compute this quantity from the correlator 
(0), we adopt a rational parametrization for Q, 



Q 



1 z 
Z 1 



1 

-1 



1 z 

Z 1 



Inserting this parametrization into the field theory ac- 
tion (||), and doing the functional integral in Gaussian 
approximation (tree level), we obtain an = a^ with 

4 = vD x , 

which is the result expected from quasiclassical transport 
theory. The weak localization correction to (rn(x, x')) 
arises from one- loop graphs of the kind shown in Fig. 1; 
see |^2| . The basic element of this graph is a 4- vertex 
representing the fourth-order term in the Taylor expan- 
sion of S with respect to Z, Z. A double line oriented 
by an arrow stands for the bare propagator (ZZ)q. All 
one-loop graphs are composed of three propagators and 
one 4- vertex. Although these graphs appear as a calcula- 
tional device for organizing the field-theoretic perturba- 
tion expansion, they do have a direct physical meaning, 
as follows. Each of the two single lines in Fig. 1 stands 
for a Feynman path contributing to the Gorkov Green 
function G R (x, x'; ±E). Double lines represent sums of 
impurity ladders with an arbitrary number of Andreev 
scattering events inserted. It is seen from Fig. 1 that one 
of the two Green function lines proceeds directly from 
the point x' to the point x, whereas the other one makes 
an excursion in the form of a double loop. The propaga- 
tor associated with the double loop is called the D-type 
cooperon What is essential here is that the charge 
of the quasiparticle during the second looping is exactly 
opposite to the charge during the first looping. This fea- 
ture makes the D-type cooperon stable with respect to 
disorder averaging irrespective of the orbital coupling to 
a magnetic field, by canceling the Aharonov-Bohm phase 
§ A-d\ accumulated in the loop. Fig. 1 also indicates the 
fact |h]] that the present variant of the weak localization 
phenomenon already affects a single Green function and 
thus the density of states. 
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FIG. 1. One-loop diagram contributing to the correlator 
(tzz(x, x')}. The electric charge of the quasiparticle during the 
second looping (—a) is opposite to the charge during the first 
looping (a). 

By evaluating the one-loop graphs in a similar manner 
as in Ref. [|l2|, we obtain 

Di f d 3 k / \ - 1 

Sau = __i Re y __ [d^ + D±k l + 2iE) . (5) 

The full spin conductivity is an = crP, + Sau + • ■ ■• Note 
that the correction is formally similar p4| to that for 
class Al, except that it explicitly depends on energy. In 
fact, it disappears with increasing excitation energy, or 
temperature, in agreement with the fact |T^ | that moving 
up in energy causes crossover from class C to class A, 
where weak localization is absent. In dimension d < 2 
the integral over wave numbers is cut off in the infrared 
by the inverse of the dephasing length L v due to inelastic 
(or quasielastic [^5| ) scattering, while for dimension d > 2 
it is UV-regularized by the inverse elastic mean free path. 

Next recall that the quasiparticle spin is assumed to 
be conserved, which allows to consider the sectors with 
spin up (s = +1/2) and spin down (s = —1/2) sepa- 
rately. Turning on the Zccman coupling is equivalent to 
shifting the excitation energy E — > E — s/iB. As a result, 
the energy dependence of the weak localization correc- 
tion translates into a field dependence. Note that this 
effect differs from weak localization in disordered metals 
PH , where the orbital coupling to a magnetic field causes 
class Al to cross over to class A. In that case, the field 
scale is set by Bo = (eDr^)^ 1 with t v being the dephas- 
ing time. In the present case the relevant field scale is 

To compute the thermal conductivity n at temperature 
T, we use the relation 

ku = V I a H (E - Sl iB) ^(E) E dE , 

S =±1/2 J ° 01 

where fr(E) = (1 + e B / T ) _1 is the Fermi-Dirac distri- 
bution, and our unit of temperature is such that fcs = 1. 
If the energy dependence of a can be neglected in the 
range < E < T, and if T < Max(/iB,r v ), i.e. 5a 
is cut off by the Zeeman energy or the dephasing rate 
r v = TjT , rather than by the temperature, we may pull 
out au(E) from under the integral sign, thus obtaining 
an analog of the Wiedemann-Franz law: 

j, = yowff) • (6) 

Here we have combined the spin up and spin down con- 
tributions, by assuming the quasiclassical term a^ to be 
unaffected by the Zeeman splitting. 



IV. ISOLATED VORTICES 

We now specialize to an extreme type-II s-wave su- 
perconductor in a weak magnetic field (but well into the 
mixed state so that the field is approximately homoge- 
neous), where quasiparticles are bound to a dilute array 
of vortex cores and the amplitude to hop from one vortex 
to another is negligibly small. In this case the problem re- 
duces to a set of decoupled one-dimensional theories, one 
for each vortex, and we formally set D± = 0. The pa- 
rameters of the one-dimensional nonlinear a model were 
calculated by solving the Usadel equation for a single 
vortex in [Q, where we found D» = C2Vp£/3Ci, and 
v J d 2 x± = 2^n7t£ 2 Ci if the integral extends over the 
area occupied by one vortex. The parameter £ is the 
dirty coherence length, ^ is the density of states of the 
normal metal, and C± = 3.16 and C'2 — 1.20 are nu- 
merical constants dependent on the vortex profile. Us- 
ing the fact that the total number of vortices equals 
the transverse area of the sample divided by half the 
square of the magnetic length Is = \J2/k /eB, we obtain 
(tJ| = 4ttC2Vn(£,/Ib) 2 vf£/3 for the quasiclassical limit of 
the spin conductivity. The weak localization correction 
is given by 

2-Dii f dk -1 
Sa \\ = j^Re —(D ll k 2 +T v + 2i(E-s t iB)) , 

where inelastic events were incorporated by shifting the 
denominator by This result applies when the de- 
phasing length L v = y/D\\/T v is shorter than the vortex 
length L|| . In the opposite, mesoscopic regime (Ln <C L v ) 
the weak localization effect was worked out in [[12| . 

The low-temperature behavior of the thermal conduc- 
tivity depends on how T v varies with T. If we as- 
sume a power law T v ~ T p with exponent p < 1 p5| , 
then an becomes constant in the energy range where 
dfT(E)/dT is essentially different from zero, and we get 
the Wiedemann-Franz law (Q) with 

a,, (jmB) = a\ - {^l%)- l Re^D^J~+i^B) . 

In the high-field regime /J.B > T v the weak localiza- 
tion correction to the thermal conductivity is cut off by 
the Zeeman splitting, giving a characteristic dependence 
Sku/T ~ -1/yB. On the other hand, if p > 1 then the 
relevant low-T regime is T 3> T v , and the weak localiza- 
tion effect is cut off by T for low fields. In that case, one 
finds K\\/T= 2^cr°- l^/f (-v/2-l)C(3/2)L T /7r/| ) i.e. the 
quantum correction is determined by the thermal length 
L T ^^/D{JT. 

The above considerations apply to a vortex array in 
the dilute limit near H&. As the field is increased, 
the quasiparticle hopping rate between vortices in an s- 
wave superconductor grows strongly. When the field is 
tuned close to H C 2 , where the system of vortex cores be- 
comes dense, the diffusion constant D± gets large and 
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the anisotropic field theory (m) three-dimensional. Since 
the quasiparticle states of the weakly disordered three- 
dimensional system are extended, a derealization tran- 
sition must take place with increasing field. Note that 
this transition is not in a new universality class, as the 
breaking of spin rotation invariance by the Zeeman cou- 
pling reduces class C to class A |l3| . Nevertheless, the 
occurrence of such a derealization transition may be of 
experimental interest, for it can be observed by varying 
the magnetic field (instead of the disorder strength or the 
chemical potential). 

V. WEAK LOCALIZATION IN THE CUPRATES 

We now adapt our results to the very interesting case 
of quasi two-dimensional d-wave superconductors such 
as the cuprates. As was stated before, the low-energy 
quasiparticles of a dirty cZ-wave superconductor in zero 
magnetic field belong to symmetry class CI. Weak local- 
ization effects in that class arise from two distinct modes 
of quantum interference || : the cooperon of type A, and 
the cooperon of type D. The former is the natural ana- 
log of the cooperon mode well known from the theory 
of disordered metals pij . When time reversal symmetry 
is broken by a magnetic field penetrating the supercon- 
ductor, the A-type cooperon becomes massive and dis- 
appears over a scale given by Bo — (eDr^)^ 1 . This 
crossover takes class CI into class C, while leaving weak 
localization due to the D-type cooperon intact. As we 
have seen, the latter mode is cut off only by the Zeeman 
energy, which becomes effective over the characteristic 
field scale £?z = (M r v) _ • Using /i = 2/jb = e/m and 
D <~ hwljm we see that the two scales are separated by 
a large factor: Bz/Bo ~ k-pi, i.e. the elimination of 
the A-type cooperon by the orbital coupling to the mag- 
netic field takes place at much smaller fields than does 
the removal of the D-type cooperon by the Zeeman en- 
ergy. This justifies our explicitly retaining the Zeeman 
coupling, while burying the orbital coupling via the intro- 
duction of a maximum entropy model. In the following, 
we take the magnetic field to be applied along the c-axis, 
and assume the system to be well into the mixed state so 
that the field is approximately homogeneous. 

The cuprates are highly anisotropic materials, consist- 
ing of weakly coupled CuC-2 planes, for which D\\ = 
D c <C D a i, = D±. At weak interlayer coupling, the 
continuum approximation leading to (^|) is not justi- 
fied in the c-direction, and we need to restore the dis- 
crete layer structure. This is done by making the re- 
placement D,,k?< — > 2t c (l — cos fey a), where a is the dis- 
tance between layers and t c is the interlayer hopping rate. 
Then, by performing the integral in (|^) over the domain 
L^ 1 < \k±\ < I and —ir/a < hi < n/a, we obtain 

Sa ± = -(2 7 r 2 a)- 1 Reln ^F s (T)/F s (T v )^j , (7) 



F s (e) = ^e + 4t c + 2i(E-snB) + y/e + 2i(E - sfiB) , 

where T = D±/£ 2 and Y v — D±/I/Z, are the elastic and 
inelastic scattering rates. 

To evaluate the consequences of this general for- 
mula, one needs to distinguish cases. For brevity, 
we concentrate on the limit defined by the condition 
that elastic scattering sets the largest energy scale: 
T > Max(4i c , 2E, fiB). Consider first the case 4t c < 
Ma,x(2E, fiB), which physically means that the coher- 
ence of the quantum interference modes is destroyed be- 
fore quasiparticles have a chance to hop between lay- 
ers. The layers then effectively decouple, yielding a two- 
dimensional system, and the formula for 5cr± becomes 

Sa x = -(47r 2 a) _1 Rcln (r/(T v + 2iE - 2isfiB)"j . 

The appearance of a logarithm is characteristic of weak 
localization in two dimensions. For the in-plane thermal 
conductivity we get 

5k ± (B)/T= -{12a)- 1 RelnTr/fT^ + iM-B)) , 

provided that the conditions of validity of the Wiede- 
mann-Franz law (^) are satisfied. Note that in contrast 
with three-dimensional metals, where weak localization is 
a rather minute effect, the correction here can easily ex- 
ceed 10% under experimental conditions. This is because 
the relative size 5k/ k is roughly given by the inverse of 
the dimensionless intralayer coupling constant 2-KV2dD\_ , 
whose value in zero field has been estimated @,^3| to be 
not much in excess of unity. 

In the opposite limit, where 4i c is much larger than 
fiB and T Vl but still smaller than T, the field depen- 
dence of the weak localization correction to the thermal 
conductivity becomes three-dimensional: 

S -^± = ~ (in ^YjTc - RcyV^ + ^B)/4t^j , 

where again the law (^) was assumed. Note that the 
above expressions for Sk±(B) increase with B. 

To summarize, weak localization in class C, cut off by 
the Zeeman splitting, causes the thermal conductivity to 
increase with the magnetic field at sufficiently low tem- 
peratures. To make this more quantitative, we need to 
specify the field/temperature range where the effect be- 
comes observable. The answer is provided by the value of 
the spin magnetic moment of the electron (with a g-factor 
of 2), which is 1.35 K/T in suitable units. As a result, 
if the field strength is of the order of 1 Tesla, the weak 
localization induced field dependence sets in at temper- 
atures below one 1 Kelvin (unless for some unexpected 
reason the dephasing rate T v is anomalously large) . 

We now wish to elucidate whether such an effect might 
already be visible in recent experiments. The discussion 
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is somewhat complicated by an ongoing debate concern- 
ing the leading, quasiclassical term k , Let us summarize 
the current situation as we see it. 



Krishana et al. 27 measured the magnetic field depen- 
dence of the thermal conductivity in a BSCCO system for 
temperatures T > 6 K. After an initial decrease at weak 
fields, they observed a sharp kink at B* ~ \/T, followed 
by a wide plateau for B > B* . The nonanalyticity at B* 
has been interpreted ]28| as a phase transition to a new 
ground state with a secondary id xy order parameter. We 
will not be particularly concerned with that issue here. 
(The addition of an id xy component to the order parame- 
ter is fully compatible with the symmetries of class C and, 
if disorder is present, the field theory (||) for the quasipar- 
ticle excitations remains qualitatively unchanged.) From 
the observation of field independence over a sizable range 
of temperatures, one deduces Q that both the electronic 
and the phonon contribution to the thermal conductivity 
must be individually constant. The constancy of the elec- 
tronic part was initially attributed to the d x 2_ y 2 + id xy 
state being fully gapped, i.e. to the complete absence 
of low-energy quasiparticle excitations. This explanation 
has been challenged by experimental data of Aubin et 
al. |29|. While confirming the results of Ref. pTl] for 
T > 5 K, these data reveal the emergence of a posi- 
tive thermal magnetoconductance at lower temperatures 
T < 1 K. (The data also show pronounced hysteresis 
effects whose interpretation remains controversial.) Tak- 
ing the constancy of the phonon contribution for granted, 
the observation of such dependence strongly indicates a 
residual density of quasiparticle states at zero energy. 
The existence of such states is no surprise. Indeed, in 
the mixed state of a superconductor with d x 2_ y 2 wave 
symmetry a residual density of states is expected even in 
the absence of disorder, because some fraction of the low- 
energy quasiparticles (close in momentum to the d-wave 
nodes) are Doppler shifted to zero energy by the super- 
current circulating around the vortices (the Volovik effect 
H), which leads to v{E = 0,B) = m\fB/B c2 . (For 
a recent discussion of the same effect for a ground state 
with d x 2_ y 2 +id xy symmetry, see The residual den- 

sity of states created by this mechanism is approximately 
constant in energy below the average Doppler shift scale 
Eb, roughly estimated by Eb/T c ~ y/B/B C 2. Disorder 
can only broaden the range of energy independence of 
v. Hence, assuming B C 2 ~ 100 T and a superconducting 
transition temperature T c ~ 100 K, the energy scale Eb 
is of order 10 K for fields of a magnitude of about 1 T. 

Now recall the experimental observations reported in 
p9| : an electronic thermal conductivity which is indepen- 
dent of the magnetic field for T > 5 K (and H > H*), and 
begins to increase with B below T ~ 1 K. (According to a 
footnote in p9[ ], the same effect has been seen in YBCO.) 
The first point to address is the field independence at the 
higher temperatures. Franz p2fl has recently proposed a 
model for the quasiclassical thermal conductivity k , in 



which the increase of v with the field is exactly canceled 
by a concomitant decrease of the quasiparticle mean free 
path I. The model assumes scattering from the super- 
flow due to randomly positioned vortices. In contrast, 
another recent theory argues in favor of the domi- 
nant scattering mechanism being impurities close to the 
unitarity limit. We will not pursue here the discussion 
as to which is the correct model to use. With the mi- 
croscopic theory of the plateau effect being a subject of 
debate, our philosophy is to accept it as an experimen- 
tal fact that the field variation of v and £ is such as to 
cancel in k° ~ v(B)l(B). The question to address, then, 
is why a field dependence sets in when the temperature 
is lowered. We argue that this is at least in part due 
to weak localization. As we have seen, weak localiza- 
tion in the mixed state of dirty d-wetve superconductors 
is a phenomenon on safe theoretical ground, is sizable in 
magnitude, and is expected to occur at the right tem- 
perature and field scales to match the experiment [p9| . 
To preclude any confusion, we stress that the effect un- 
der consideration is distinct from weak localization in 
combination with Aslamasov-Larkin fluctuations, which 
have been invoked in |34j to explain the negative thermal 
magnetoconductance observed in a dirty LSCO system at 
much higher temperatures. 

Theories proposed by previous authors attribute the 
temperature variation of dn(B,T)/dB to the leading 
(quasiclassical) term, n° ~ vt. Given the low-energy 
constancy of the density of states, such a variation would 
have to arise from an energy (or temperature) depen- 
dence of the elastic mean free path. Possible explana- 
tions are: i) low-energy transparency of d-wave vortices 
to quasiparticles p2] |, and ii) energy-dependence of the 
elastic scattering rate due to impurities near the unitarity 
limit ^33|. The challenge to these scenarios is to explain 
why for fields B ~ 1 T the effect sets in at tempera- 
tures around 1 K. In the weak localization scenario we 
have described, this comes about very naturally if T v is 
determined by thermal broadening, since fi = 1.35 K/T. 

A clear difference is that weak localization effects con- 
tinue to be enhanced with decreasing temperature - they 
ultimately drive the system to an insulator by localiza- 
tion of all quasiparticle states - whereas the energy de- 
pendence of the elastic mean free path saturates. To 
discriminate, it is therefore desirable to push the exper- 
imental measurements to the lowest temperatures possi- 
ble. In order to achieve a quantitative description based 
on formula (Q), it will be necessary to take the field de- 
pendence of £ into account. Our suggestion is to extract 
the density of states v(B) from measurements of the spe- 
cific heat, and then deduce £(B) from the quasiclassical 
formula n(B) ~ u(B)£(B), valid at high temperatures 
(T > 5 K). As far as the temperature dependence of 
the dephasing rate r^, is concerned, a phenomenological 
model needs to be used. To our knowledge, a theory for 
this quantity in the mixed state of dirty d- wave supercon- 
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ductors does not exist. In the long run, weak localization 
may turn out to be the appropriate tool to measure T v , 
as is established practice in disordered metals fl24,Ea] . 



VI. CONCLUSION 

Noninteracting electrons subject to disorder and a 
magnetic field are well known to belong to the standard 
universality class A (unitary symmetry, (3 — 2). When 
spin-singlet pairing correlations are added, the univer- 
sality class of the low-energy quasiparticles changes to 
type C . It has been shown that the transport proper- 
ties of these quasiparticles are unconventional. In par- 
ticular, there exist modes of destructive quantum inter- 
ference which survive the orbital coupling to a magnetic 
field. They are cut off at higher fields by the Zeeman cou- 
pling, thereby giving rise to a field dependent quantum 
(or weak localization) correction to the low temperature 
thermal conductivity, with the characteristic scale given 
by [i = 1.35 K/T. A good place to look for such correc- 
tions experimentally are disordered low-dimensional su- 
perconductors, such as the cuprates, in the mixed state. 

On general symmetry grounds, the low-energy effective 
field theory for quasiparticles in class C is predicted to be 
a nonlinear a model of type £)III|CI. The Lagrangian of 
this field theory has a universal form, independent of the 
symmetry of the order parameter (s, d, etc.), as long as 
the superconductor conserves the quasiparticle spin and 
is penetrated by magnetic flux. The role of the supercon- 
ducting ground state is merely to determine the values of 
the field theory coupling constants, their anisotropy, and 
their dependence on energy and magnetic field. Quanti- 
tative predictions for the weak localization corrections to 
transport can be made once the values of the couplings 
and their dependences have been obtained, either from 
quasiclassical transport theory or from experiment. We 
advocate the use of such predictions in understanding the 
low-temperature experiments of Aubin et al. 
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